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Abstract 

The main purpose of these lectures is to discuss briefly recent meth- 
ods of calculation of statistical properties of quantum eigenvalues for 
chaotic systems based on semi-classical trace formulas. Under the 
assumption that periodic orbit actions are non-commensurable it is 
demonstrated by a few different methods that the spectral statistics 
of chaotic systems without time-reversal invariance in the universal 
limit agrees with statistics of the Gaussian Unitary Ensemble of ran- 
dom matrices. The methods used permit to obtain not only the lim- 
iting statistics but also the way the spectral statistics of dynamical 
systems tends to the universal limit. The statistics of the Riemann 
zeta function zeros is considered in details. 
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1 Generalities 



1.1 Trace formulas 

The quantum density of states for a chaotic dynamical systems 

d{E) =Y^5{E - E„) (1) 

n 

under quite general conditions can be expressed through a sum over classical 
periodic orbits by the Gutzwiller trace formula (plus corrections 
if necessary) 

= + c/(°^")(E), (2) 

where the smooth part, d{E), for a /-dimensional system with a Hamiltonian 
H{p,q) is given by the Thomas- Fermi formula 

and the (most interesting) oscillating part of the level density, d^°'^'^\E), has 
the following form 

1 oo 1 c 

rf(°-)(^) = y T„y — -— -_cos(n(^ - -/i„)). (4) 

TT^i^ ^^Jdet(Mp"-l)|V2 V ^^-P)) KJ 

Here the first summation is done over primitive classical periodic orbits (ppo) 
with enrgy E and the second one is performed over all repetitions of a given 
ppo. Sp is the classical action calculated for a ppo labeled by p, Mp is the 
monodromy matrix for this trajectory, and Hp is its Maslov index [Q. 

Similar formulas can be written also for other spectral functions. 
In particular the staircase function has the following form 

N{E) = d{E')dE' = N{E) + Ar(°^'=)(E), (5) 

J —oo 



where 



{2'nK) 



and 

1 oo -1 C _. 



^;ron- ^|det(M«-i)|i/2 ^ 2' 
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We shall need also the dynamical zeta function which (for 2- dimensional 
systems) equals the following product over ppo 



oo piSji/'h—i-nfip/2 

Here Ap is the largest eigenvalue of the monodromy matrix (|Ap| > 1) 

MpU = ApU. (9) 

This zeta function is connected with the density of states and the staircase 
function by the following relations 

rf(-)(i^) = _^(^_c.c.), (10) 

and 

Ar(°^'=)(E) = -^(InZ(E) - c.c). (11) 

The trace formulas exist not only for dynamical systems but also for the 
Riemann zeta function (and others number-theoretical zeta functions 
as well). 

The Riemann zeta function is defined as 

a^)-t^.-^J^y (12) 



n=l 



where the product is taken over prime numbers. This function converges only 
when Res > 1 but it is well known |^ that it can be analytically continued 
in the whole complex s-plane, the only singularity being the pole at s = 1 
with unit residue. 

According to the famous Riemann conjecture all nontrivial zeros of this 
function have the form Sn = \ + iEn and the density of these zeros d{E) = 
J2n ~ En) can be expressed by the following 'trace' formula 

d{E) = d{E) + S"''\E), (13) 

where 

JW^^ta^. (14) 
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and 

1 °° In n 

S-\E) = --EE 3|cos(Enlnp), (15) 

p n=l " 

where the summation is performed over all prime numbers. 

The staircase function for the Riemann zeros has a similar form 

N{E) = N{E) + N'-"'''\E), (16) 

where 

^W = |(>n|-l)H-5. (17) 

and 

1 oo 1 

^(o.c).^N yy ^sin(Enlnp). (18) 

By comparing Eq. (^ and Eq. (|I^) one observes a remarkable correspondence 
between different quantities in these trace formulas 

Periodic orbits of chaotic systems primes, 
periodic orbit period In p. 

The convergence properties of both formulas are also quite similar. The 
number of periodic orbits with period less than T for chaotic systems is 
asymptotically (see e.g. 0) 

hT 

N{T, <T) = —, (19) 

where constant h is called the topological entropy. The number of prime 
numbers less than x is given by the prime number theorem |^ 

N(p <x)= n(x) = (20) 
mx 

As Inp = Tp this expression has the form similar to (|19D with h = 1 

N{Tp < T) = ^. (21) 

Due to these similarities number-theoretical zeta functions play the role of 
simple (but non-trivial) models of quantum chaos and in these lectures we 
consider in parallel both dynamical systems and zeta functions and shall 
check our methods first on the Riemann case. 
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1.2 Random matrix theory 

Wigner and Dyson in the fifties liad proposed to describe complicated (and 
mostly unknown) Hamiltonian of heavy nuclei by a member of an ensemble 
of random matrices (see |0) and they argued that the type of this ensemble 
depends only on the symmetry of the Hamiltonian. 

For systems without time-reversal invariance the relevant ensemble is the 
Gaussian Unitary Ensemble (GUE), for systems invariant with respect to 
time-reversal the ensemble should be the Gaussian Orthogonal Ensemble 
(GOE) and for systems with time-reversal invariance but with half-integer 
spin energy levels should be described according to the Gaussian Symplectic 
Ensemble (GSE) of random matrices 0. For these classical ensembles 
all correlation functions can be written explicitly [0. In particular 2-point 
correlation functions 

R2{x) = l + 5{x)-Y2{x), (22) 



have the following form [0, [T^ 
For GOE 



y-A.^) = (—? - (S.(-) - - 7^)- (23) 



Here 



and e{x) = sgn(a;). 
For GUE 



For GSE 



Si(a:) = £ '^dy, (24) 



Y2{x) = ^J^y. (25) 
nx 



, ,sin27rx,n ^^cos27rx sin27rx, 

2-point correlation form factor is the Fourier transform of the 2-point 
correlation function ^ 

K{t) = [ R2{x)e^^''''dx, (27) 



and for three classical ensembles it has the following forms. 
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For GOE 



r 2t-tln(l + 2t) 0<t<l . . 

1 2-tln((2t + l)/(2t- 1)) t>l ■ 



For GUE 




(29) 



For GSE 

A-W = | Y-'"—'" . (30) 

The nearest neighbor distributions for classical ensembles can be expressed 
through solutions of certain integral equations and numerically they are close 
to the Wigner surmises 

GOE: p(s) = ^.exp(-p). (31) 

GUE: p{s) = %s\xp{ s^). (32) 

GSE: p(^) = T^^'exp(--.^). (33) 

Later it was understood that the same conjectures can be applied not only 
for heavy nuclei but also for simple dynamical systems and to-day standard 
accepted conjectures are the following 

• The energy levels of classically integrable systems on the scale of the 
mean level density behave as independent random variables and their 
distribution is close to the Poisson distribution (Berry, Tabor P]). 

• The energy levels of classically chaotic systems are not independent but 
on the scale of the mean level density they are distributed as eigenvalues 
of random matrix ensembles depending only on symmetry properties 
of the system considered (Bohigas, Giannoni, Schmit 0). 

— For systems without time-reversal invariance the distribution of 
energy levels should be close to the distribution of the Gaussian 
Unitary Ensemble (GUE) characterized by quadratic level repul- 
sion. 
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— For systems with time-reversal invariance the corresponding dis- 
tribution should be close to that of the Gaussian Orthogonal En- 
semble (GOE) with linear level repulsion. 

— And for systems with time-reversal invariance but with half-integer 
spin energy levels should be described according to the Gaussian 
Symplectic Ensemble (GSE) of random matrices with quartic level 
repulsion. 

These conjectures are very well confirmed by numerical calculations (see e.g. 
1^, fill). The main purpose of these lectures is the developing of methods 



which permit to prove analytically one part of these conjectures, namely that 
in the universal limit (i.e. on the scale of the mean level density) the 2-point 
correlation function for systems without time-reversal invariance agrees with 
Eq. (|25D. For further references it is convenient to rewrite this expression in 
the dimensional form 

^2(e) = (Pi?2(rfe). (34) 

It gives 

R,{e) = + R^ie) + Rt'\e), (35) 



where the smooth part 
and the oscillating part 



R2{e) = (36) 



p2mdt _|_ p—2mdt 

R'r'(^) = ■ (37) 

The plan of the paper is the following. In Section |^ we relate through the 
trace formulas correlation functions of quantum eigenvalues (and of Riemann 



zeros) with periodic orbit sums and in Section |2.1| the simplest approxima- 



tion of evaluating these sums called the diagonal approximation is discussed. 



But as shown in Section |2.2| this approximation is valid only for large energy 
difference and to find correlation functions in the full range new methods 
are needed. In Section |^ a method specific for the Riemann zeta function is 
discussed. To calculate off-diagonal terms we use the Hardy-Littlewood con- 
jecture for distribution of pairs of primes and in Section |0] it is demonstrated 



that this conjecture leads to interesting formulas for the 2-point correlation 



7 



function of the Riemann zeros. In the universal hmit this correlation func- 
tion tends to the GUE results but it also permits to investigate how the 
correlation function tends to the universal result. In Section P]l we briefly 
mention a certain model where a generalization of the Hardy-Littlewood con- 
jecture also permits to calculate the 2-point correlation function. But the 
Hardy-Littlewood conjecture can be aplied only for primes and it cannot be 
generalized for dynamical systems. In Section ^ a certain method is pro- 
posed which overcomes this difficulty. The main ingredient of this method 
is an artificial construction of approximate density of states which has the 
correct analytical properties but requires the knowledge only of finite num- 
ber of periodic orbits. In Section ^TT] it is shown that for the case of the 
Riemann zeros the method gives exactly the same result as was obtained by 
using the Hardy-Littlewood conjecture. In Section ^]2|this method is applied 
to dynamical systems and it was demonstrated that under the assumption 
that all periodic orbit actions are non- commensurable it is possible to prove 
that spectral statistics of generic dynamical systems without time-reversal 
invariance in the universal limit tends to the GUE statistics. In Section ^ an 
other method of calculation of spectral statistics is proposed. The method 
is based on the universality of different random matrix ensembles and it is 
shown that it gives exactly the same expressions for 2-point correlation func- 
tions for the Riemann zeros and dynamical systems as have been obtained 
in previous Sections by different methods. Finally in Section ^ one more 
method is discussed which is based on the Riemann-Siegel type resummation 
of the trace formulas which also leads to the same expressions for 2-point 
correlation functions. Though the exact mathematical proof of the above re- 
sults is not known and all our methods should be considered as heuristics or 
'physical' proofs the agreement between different methods strongly indicates 
the correctness of the results. 



2 Correlation functions 

Formally ra-point correlation functions of energy levels are defined as the 
probability of having n energy levels at given positions and they are connected 
to the density of states by the relations 

K(ei, £2, . . . , en) =< d{E + ei)d{E + £2) . . . d{E + e„) >, (38) 
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where the brackets < . . . > denote the smoothing over an appropriate energy 
window 

< fiE) >= J f{E')a{E - E')dE', (39) 

with 

j a{E)dE = 1. (40) 

Function (j{E) is assumed centered near zero and has a width /\E which 
obeys inequahties 

AEq < < AEe/ < E. (41) 

Here AEq is of the order of the mean level separation, /S.Eq 1/d, and IS.Eci 
denotes the energy scale at which classical dynamics changes noticeable. 
Let us write the trace formula for the density of states (^ in the form 

d{E) = d{E) + J2 TpAp,„e'"^''(^)/^ + c.c, (42) 

p,n 

where 

27r^|det(M«-l)|V2^ " ^""^^ 

Substituting this expression to the formula for 2-point correlation function 
one gets 



1 

Yl ^Pi^P2^pi,ni^p2,„2 < exp -{Sp,{E + ei) - Sj,^{E + eg)) > +c.c., 



i?2(ei,e2) (44) 

i 

Pl,P2 

and the terms with the sum of actions are assumed to be washed out by the 
smoothing procedure. 

Expanding the actions and taking into account that dS{E)/dE = T(E) 
where T{E) is the classical period of motion one finds 

i?2(ei,e2) = +Y1 ^Pi,ni^P2,n2 < exp-(ni5p,(E) - n2Sp^{E)) > 

X exp -{niTp^{E)ei - n2Tp^{E)e2) + c.c. (45) 
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2.1 Diagonal approximation 

Berry proposed to estimate this sum by taking into account terms only 
with exactly the same actions having in mind that terms with different values 
of actions should be small after the smoothing. Therefore in this approxima- 
tion (called the diagonal approximation) 



Rf^^\e) = Y: T^lA.^Xe'''^''^^^' + c.c, (46) 

p,n 



Here e = ei — €2 and the sum is taken over all periodic orbits with exactly 
the same action. 

Introducing the classical zeta function [] 

sT 

Zci{s) = m^-hr\^~' (47) 

ppo M^pl 

one can rewrite the 2-point correlation function in the form 



Ri'"''\^) = --;^^^nAie), (48) 



where 

A(e) = |Z,K^e)|2$(^-^')(e), (49) 

and function <|)('**'*5)(e) is defined as the following convergent sum over peri- 
odic orbits 

n) - exp (e E ^(;^^(^-L^ - l)e-.<-) + «.) . (50) 

Another useful form of these relations is the expression for the diagonal 
approximation of 2-point form factor defined in Eq. ( pTj) 

Xidiag)^^^ = 27: J2T^\Ap,n\^ 6 {27ct - nTp{E)) + c.c. (51) 

p,n 



"'^ Often the classical zeta function is introduced in such a manner that 
EppoEZi fc|dct(Mp^-i)| = In^c/, which gives Z,,(s) = UppoU^^oi^ " JK^)~^""^^^- 
As the convergent properties of both zeta functions are the same we chose the simplest 
definition sufficient for our purposes. 
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According to the Hannay-Ozorio de Almeida sum rule for ergodic systems 



12] one can compute such sums by substituting the local density of periodic 



orbits 



(52) 



and consequently 

K^'''"'^ it) = ^j Tp5{2^t - Tp)dTp = gt, (53) 

where g is the mean multiplicity of periodic orbits (i.e. the number of pe- 
riodic orbits with exactly the same action). For generic systems without 
time-reversal invariance there is no reasons for equality of actions for differ- 
ent periodic orbits and g = 1 but for systems with time-reversal invariance 
each orbits can be passed in two directions therefore in general for such sys- 
tems g = 2. Comparing ( |53D with Eqs. ( pS] ) and (^) one concludes that 
the diagonal approximation reproduces the correct small-t behavior of form- 
factors of classical ensembles. 

Unfortunately, K'^'^'^a) [t) grows with increasing of t but the exact form- 
factor for systems without spectral degeneracy should tends to d for large t 
which reflects the existence of the delta function in -R2(e) for small e [^] 

i?2(e) ^ d5{e), when e ^ 0, (54) 

or 

K{t) — * d, when t — > oo. (55) 

This evident contradiction clearly indicates that the diagonal approximation 
cannot be correct for all values of t and more complicated tools are needed 
to obtain the full form-factor. 



2.2 Criterion of applicability of diagonal approxima- 
tion 

One can give a (pessimistic) estimate till what time the diagonal approxima- 
tion can be correct by the following method. The main ingredient of the di- 
agonal approximation is the assumption that after smoothing all off-diagonal 
terms give negligible contribution. This condition is almost the same as the 
condition of the absence of quantum interference. But it is known that the 
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quantum interference is not important for times smaller than the Ehrenfest 
time which is of the order of 



tE^jH^/h), (56) 

where A is a (classical) constant of the order of the Lyapunov exponent 
defined in such a way that the mean splitting of two nearby trajectories for 
the time t grows as exp(Ai). For billiards A = kX' where k is the momentum 
and A' determines the deviation of two trajectories with the length L — kt, 
and the role of H plays k~^ 

tE^^Hk)- (57) 

More refined estimates can be done as follows. The off-diagonal terms can 
be neglected if 

< exp ^{Sp,{E) - Sp,{E)) X 1. (58) 

But this quantity is small provided the difference of period of two orbits 
AT = Tpi — Tp^ times the energy window AE used in the definition of 
smoothing procedure is large 

i(Tpi-Tj,JAE»l. (59) 

For billiards Tp = Lp/k and this condition means that one has to consider 
all periodic orbits such that their difference of lengths is 

hk 

Lp, - Lp, > (60) 
But the number of periodic orbits with the length L grows exponentially 

N{Lp <L)^—, (61) 

where h is of the order of the Lyapunov exponent A'. Therefore in the interval 
L,L + 61 there is e^^5l/L orbits and the mean difference of lengths between 
orbits with the length < L is of the order of 

AL = Lexp{-hL). (62) 
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To fulfilled the above condition one has to restrict the maximum length of 
periodic orbits, Lm, by 

kh 

Lmexp{-hL„,) ^ (63) 



In the limit of large it gives 



1 , AE 



which corresponds to the Ehrenfest estimate above. 

Let us denote h = 1/Iq where Iq has the dimensionality of the length. 
Then the above estimate can be transform into the following form 

AE 

L^^ lo'^n-—, (65) 
and Et is an analog of the Thouless energy 

Et = — and tt = -r- (66) 

tt k 

Note that the Heisenberg time is 

tH = 2TTd, (67) 

and for billiards because J is a constant 

tE < tH. (68) 

For the Riemann zeta function the situation is better because the role of 
'energy' in this case plays the 'momentum' and the density of states for the 
Riemann zeta function is (ln(ii^/27r))/(27r). As in this case h = 1 

tE = tH, (69) 
and the diagonal approximation for the Riemann zeta function is valid till 



the Heisenberg time [|13|, pA 



This type of estimates directly indicates that the diagonal approxima- 
tion for dynamical systems can not, strictly speaking, be used to obtain 
an information about the form-factor for large value of t. (Note that for 
GUE systems the diagonal approximation gives the expected answer till the 
Heisenberg time but it just signifies that one has to find special reasons why 
all other terms cancel.) 
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3 Beyond the diagonal approximation 

The simplest and the most natural way of semi-classical computation of 2- 

point correlation functions is to find a method of calculating off-diagonal 
terms. Wc shall discuss here this type of computation on the example of the 
Riemann zeta function where much more information than for dynamical 
systems is available. 

The trace formula for the Riemann zeta function may be rewritten in the 
form 

^(o.c)(^) ^Hn) cos{Elnn), (70) 

where 



^ n=l 



^^^^ I 0, otherwise ' ^'^^^ 

The connected 2-point correlation function of the Riemann zeros, R2'^ — 
i?2 — ^, is 

The diagonal approximation corresponds to taking into account terms with 

Hi = 712 



n 

,2^ 



This expression may be transform as follows 

i?f-^(e) = -^^lnA(e), (74) 

where 

A(e) = |((l + ie)|2$(<i^«9)(e), (75) 
and function (^^^ jg given by a convergent sum over prime numbers 

°° 1—77? 

^{dia9),A ^ exp(V V L^Jlle^rnlnpe ^ y /yg) 
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In the limit e ^ C(l + ie) (ie) ^ and $('^*'^5')(e) — > const. Therefore in 
this hmit 

(77) 



which agrees with the smooth part of the GUE resuh (|36D. 
The off-diagonal contribution has the form 

R^°^-'^\ei €2) = — y ^(^1)^(^2) ^ ^ ^»£:in(ni/n2)+»(6i Inm-ea Inna) ^ ^^^^ 

(78) 

The term exp{iEln{ni/n2)) oscillates quickly if ni is not close to n2- Denot- 
ing 

rii = n2 + d (79) 
and expanding all smooth functions on d one gets 

R'2"\e) = ^ E ^-i-^-^^- > +C.C.), (80) 

where e = ei — €2- 

The main problem is clearly seen here. The function 

F{n,d) = A{n)A{n + d) (81) 

is quite a wild function as it is nonzero only when both n and n + d are power 
of prime numbers. As we have assumed that d <^ n it is quite natural to 
assume that the dominant contribution to the 2-point correlation function 
will come from the mean value of this function over all n, i.e. one has to 
substitute into R2^^\e) instead of F{n,d) its mean value 

1 ^ 

a{d)= lim — V A(n)A(n + rf). (82) 

n=l 

3.1 The Hardy-Littlewood conjecture 

Fortunately the explicit expression for this function comes from the famous 
Hardy-Littlewood conjecture ||15[. There are two different forms of this con- 
jecture which, of course, can be mutually transformed 

p — 1 



a 



id) = C2ll^^ (83) 



p\dP 
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where the product is taken over all prime divisors of d bigger than 2 and C2 
is the so-called twin prime constant 

C2 = 2 n(l-^^TT^)~ 1-32032.... (84) 

The other form which will be useful for us is expressed through the so-called 
singular series 

where the sum is taken over all natural q and all integer p co-prime to q {p < 
q). Function /x(n) is the Mobius function defined through the factorization 
of n on prime factors 

1 if n = 1 

= ^ (-1)^ ifn = pi...pk . (86) 

if n is divisible on 

Function ip{n) is the Euler function which counts the number of integers 
smaller than n and co-prime to n 

V'(n)=nn(l--), (87) 

p\n P 

where the product is taken over all prime divisors of n. 
Taking the above formulae as granted we get 

R2"\e) ^ ^j:-e'^'--j:a{d)e'^i +C.C. (88) 

After substitution the formula for a{d) and performing the sum over all d 
one obtains 

R^2"\e) = -L^y:-e-^- E (^Xs{^-^) + c.c., (89) 

where the summation is taken over all pairs of mutually co-prime positive 
integers p and q (without the restriction p < q). 
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Changing the sum over n to the integral permits to transform this ex- 
pression to the sum over that values of n where 

q /vrn 

and in this approximation 

Rr\e) = -i,e--^ E + cc. (90) 

Using the formula 

oo 

E /(p) = EE/(^'^)M^)> (91) 

ip,q) = l k=l S\q 

and taking into account that 271(1 = ln(i?/27r) we obtain 

4°^^)(e) = -i-|C(l + te)\^e^''"''^^°f^\e) + cc, (92) 
where function (e) is given by a convergent product over primes 

$(°//)(e) = n(l-^^^), (93) 

and <l>(°^^)(0) = 1. 

In the limit of small e 

= 777^(6''^^'^ + e-'"''"^), (94) 

which exactly corresponds to the GUE results for the oscillating part of the 
2-point correlation function (p7|). 

The above calculations demonstrates how one can compute the 2-point 
correlation function through the knowledge of pair-correlation function of pe- 
riodic orbits. For the Riemann case one can prove under the same conjectures 
that all ra-point correlation functions of Riemann zeros tend to corresponding 
GUE results 0. 
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The interesting consequence of the above formula is the expression for 
the 2-point form-factor 



This formula means that the off-diagonal 2-point form factor is a sum over S- 
functions in special points which are situated in a vicinity of the Heisenberg 
time plus a difference of periods of two pseudo-orbits (= logarithm of the 
difference between two integers). This set of 5-functions is dense but the 
largest peaks correspond to the shortest pseudo-orbits. Similarly the 2-point 
diagonal form factor is the sum of 6 functions in the positions of periodic 
orbits 

A--'W^jiji:!^^(*-^lnp). (96) 

The smooth values corresponding to the random matrix predictions appears 
only after a smoothing of these functions over a suitable interval of t. 



3.2 Arithmetical systems 



Similar behavior has been observed in a completely different model, namely 
for distribution of eigenvalues of the Laplace-Beltrami operator for the mod- 
ular domain |T^. Using a generalization of the Hardy-Littlewood method 
it was shown that in this model the 2-point correlation form factor can be 
written in the following form 



— E 



p 



(97) 



where 



and 



t =hn^ 
K irp 



Pip, q) 



S{p,P] q) 



(98) 
(99) 
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Here the product is taken over all prime divisors of q and S{p,p;q) is the 
Kloosterman sums 



S{n,m;c)= ^ exp ( — [nd + md (100) 

{d,c)=l ^ ^ ^ 

This model belongs to the so-called arithmetical models which are models on 
the constant curvature surfaces generated by discrete arithmetic groups. For 
all these models due to the exponential multiplicity of periodic orbits one 



expects that the spectral statistics will tend to the Poisson distribution 
though from classical point of view all these models are the best known 
examples of classically chaotic motion. Using the above expression one can 
prove this statement for the modular domain. 



4 Construction of the density of states from 
finite number of periodic orbits 

The main difficulty in using the trace formulas is their divergent character. 
They cannot converge on the real they have to produce (5-functions 

singularities there. Usually they are defined by an analytical continuation 
from a region in the complex plane of energy which do not touch the real axis. 
E.g. the Riemann zeta function converges only when Re s > 1 but zeros are 
assumed to lie on the axis Re s = 1/2. The diagonal approximation consists 
in some sense on the computing the density of states from a sum over a 
finite number of periodic orbits but such sums can never have 5-function 



singularities. In this section we shall discuss a special method |T9[ which 
permits to avoid this difficulty and produces an (artificial) expression for 
the density of states with required singularities from the knowledge of finite 
number of periodic orbits. 

Let us write the semi-classical formula not for the density of states but 
for the staircase function 

Nt,{E) = N{E) + N^T\e), (101) 

where the oscillating part of this function is truncated (smoothly if necessary) 
so to include periodic orbits with period up to a fixed period T*, where T* 
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is a parameter to be fixed later 

oo 

► \ A Rin (n ( 

n 2 



oo _ c 

N^T\e) = 2 Y: E4nsin(n(-^ --/.,)), (102) 



and for dynamical systems the pre-factor is 

1 



27r|det(M^"-l)|V2' ^^^^^ 
and for the Riemann zeta function 

In the limit T* oo function Nt*{E) should have a unit jump each time 
when E equals an eigenvalue. For finite value of T* one can at the best have 
a smooth increase at these points. The idea of the proposed method is to 
define semi-classical eigenvalues En according to the following 'quantization 
condition' |]21[] 

NT*{En)=n + ^. (105) 

The main advantage of this method is that it cannot miss any one level 
simply because lines iV = n + 1/2 will in any case cross the curve N{E) and 
will produce (approximate) semi-classical energy levels. In principle, one can 
obtain additional levels if the curve N^^iE) has decreasing parts. But for 
not too big number of included orbits it is not the case and this method is 



quite efficient for numerical computations of semi-classical levels [20 



The other important point (which explain why we take into account the 
infinite sum over repetitions of periodic orbit) is that 



2niN ^T" {E) 



exp(2vrziV^.(E)) = e'^^^' (106) 
where 

oo piSp/h—i-Kiip/2 

^T<E)= n n (1 - ^^nTix^), (107) 

is a truncated product over periodic orbits and the quantization condition 
(|105|) is equivalent to calculation of semi-classical energy levels En from the 
condition 

ZT*iEn) = 0, (108) 
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where Zt*{E) is a special form of the dynamical zeta function 

Zt-^E) = zt'{E) + e2"^^4.(E). (109) 

The above expression has a Riemann-Siegel form pH] and automatically 
obeys the important functional equation for the zeta function 

Z{E) = e^^'^Z*{E), (110) 

which is crucial in proving the correct analytical properties of the zeta func- 
tion. 

Let us define a new bootstrapped density 

DT4E) = J2^iE-En), (111) 

n 

where instead of exact energy levels we use the semi-classical ones defined as 
solutions of quantization condition ( |105| ) (or ( |108| ). Rewriting Dt*{E) in the 
form 

DT*iE) = dr* (E) (Nt* (E) -n~ 1/2), (112) 

n 

where (It* (E) = dNx* (E) / dE and using the Poisson summation formula one 
gets 

oo 

Dr-iE) = dr^iE) ^ (-l)^exp(27rifciVT*(^))- (113) 

fc=— oo 

We called it a bootstrapped density as it has all Fourier harmonics and not 
only the ones with period T < T*, the role of the effective orbits introduced 
beyond T* being to generate the correct analytical properties associated with 
the discreteness of quantum spectrum. 

Substituting this expression to the formula for the 2-point correlation 
function one gets 

i?2(ei,e2) = <dT4E + e,)dT*{E + e2)Y.i-'^)'"''" 

ki,k2 

X exp{27Ti{kiNT4E + ei) - Nt^E + ea))) > . (114) 
Consider first the fci = ^2 = term, which we write in the form 

<dT4E + e,)dT*{E + e2) >= rf^ + i^f^^^l^i, ^2), (115) 
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where 
and 

lnA(ei,e2) =47r2 < N^"r\E + e,)N^T\E + €2)) > . (117) 

Because only a finite number of orbits enters this expression we assume that 
the resulting sum over pairs of periodic orbits can be replaced by the diagonal 
terms, for which only orbits with exactly the same actions will contribute. 
The result can be expressed in the form 

A(el,e2) = |Z,(^e)|^ (118) 

where 

ZM= n U^M—9p\Ap,n\'e-^n, (119) 

Tp<T* n=l ^ 

and Qp is the number of orbits with period Tp (multiplicity of periodic orbit) . 

In generic systems the multiplicity Qp is the same for almost all orbits and 
so, if its value is denoted g 

Z,{s) = Z^(s), (120) 

where Z{s) is defined by the same formula but with Qp = 1. For systems 
without time- reversal invariance g = I, and for systems whose dynamics is 
time-reversal symmetric g — 2. 

Since the ki — k2 — term corresponds to the usual diagonal approxima- 
tion, the other terms are representing the off-diagonal contributions, R^2^^\ 
To evaluate them we note that the Taylor expansion in powers of e of the 
mean staircase function leads to a term {ki — k2)N{E) in the phase which is 
of the order of h~-^ . Hence the energy average renders negligible any contri- 
butions with ki ^ k2, and therefore 

^^^"^('^ = TTTT E 7^ eM^^idk(e, - e2))Mei, 62), (121) 
oeiOe2 ^_^Q [ZnK) 

where 

*fe(ei, £2) =< exp(27r^A;(iV^'/'^(E + ei) - N^T\e + £2))) > • (122) 
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We now make the key assumption that in generic systems the orbits up to 
period T* are non-commensurable modulo exact degeneracies, and hence that 
the energy average of any smooth function of exp{iSp) 



< />=< /(e^^i(^),e'^^(^'),...,e^- 



'Sm{E) 



(123) 



can be calculated using 




M 



d(j)j 
27 



(124) 



3 



This is essentially equivalent to random phase approximation, or to ergodic 
theorem for quasi-periodic functions with non-commensurate periods, or to 
strict diagonal approximation. 

Using the trace formula for the staircase function one can thus perform 
the energy average by evaluating these integrals. Exact results for certain 
cases will be discussed later. First, for clarity, we consider a simple leading 
order approximation to based on the relation 



This is an identity if G is a Gaussian random function with zero mean. But if 
one ignores all terms with repetitions of the same periodic orbits the formula 
for A^(°^^)(^) will be equal to a sum of a big number of terms with non- 
commensurable frequencies which can be considered as independent random 
variables therefore the resulting distribution of N'^°^'^\E) being a sum over 
many independent random variables should a Gaussian distributed random 
function. Consequently this type of approximation is expected to be a good 
approximation for generic systems. 
From it one obtains 



^fe(ei,e2) = <exp(-27rV<(4°r\£; + ei)-4°r)(£; + e2))')> 




(125) 




(126) 



where A(e) was defined above and 




(127) 
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Using the Hannay-Ozorio de Almeida sum rule 



, it may be shown that 



L ^ {ry. (128) 

Since we anticipate taking T* ^ Th, it then follows that the terms in the k 
sum decrease rapidly as (c?)"^^^ . Hence in the leading semi-classical order 
as ed — >■ cx), we may retain just the k = ±1 contributions and when deriving 
over e take into account only terms with exp{2TTide) . Finally one gets 

iJ<°"'(.,...) = <Fe-«-- < + + ^ ^^^^^ 

Zt-{E + ei)z:^,{E + 62) 

The last step consists in performing the exact average in this formula. We 
shall consider the case of the Riemann zeta function first. 



4.1 OfF-diagonal terms for the Riemann zeta function 

For the Riemann zeta function the the ratio of four truncated zeta functions 
in the above formula has the form 

< £i4l±iii£i4l±4 >=< n (130) 



where 



/ . -J. ■ {2)^ 



andoe[B 



(1 — A p^vp+^'^v Vl — A p~ 

^^('^-^ = 7 .a ' : (131) 

(1 - Ape-*<^f-*^p )(1 - Ape^-^f+^^p ) 

Ap=^, 0p = Elnp, r« = eanp. (132) 
VP 

As the logarithms of prime numbers are non- commensurable, 6*"^^ act as 
independent random variables for each p and 

< n ^p>= n <^p>- (133) 

lnp<T* lnp<T* 

The mean value of an individual Rp is 

2-K 



rZ-K 

< Rp> = / i?p(0) 





^^^'27r 

1 / (1 - Apze'^^'^ ) (1 - ApZ-^e-'^^^^ ) dz 



2m J {I - ApZ-^e-'^p'''){l - Apze'^p''^) ^ 
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(134) 



where the integral is taken over a unit circle in the complex z-plane. 

The last integral 
inside the contour. T' 
Simple algebra gives 



The last integral is easily computed by the residues. There is two poles 

_■ (1) 

inside the contour. The first is at 2; = and the second one is dbt z — -^7>6 



where Xp = r^^) — . 

The total contribution equals the product over all primes up to T* 

i?2(e) = C2exp(27riJe)|C(l + ze)|2$(°//)(e), (136) 

where 
and 

C = dX{{l-\). (138) 

All products in these expressions include prime numbers up to Inp = T*. 
The two first products converge when T* ^ 00 and only the last one require 
regularisation. But our parameter T* has not yet been fixed. Let us choose 
it in such a way that 

It is easy to check [Q] that asymptotically this equation leads 

T* = 27rJe^, (140) 

where 7 is the Euler constant. Therefore, our T* is of the order of the 
Heisenberg time as was expected. Note that exactly the same factor e'^ often 
appears in the statistical approach to prime numbers (see discussion in |T3|] ) 
and can be considered as a renormalisation of formally divergent sums. 

After this renormalisation we get exactly the same formula as has been 
derived in the previous section using the Hardy-Littlewood conjecture about 
the pairwise distribution of prime numbers. Note that in a present derivation 
no analogous conjectures have been assumed. 
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4.2 OfF-diagonal contribution for dynamical systems 

Let us now compute off-diagonal contribution to dynamical systems. Our 
starting point will be the same expression as for the Riemann zeta function 

The only difference with Riemann case is that the truncated zeta function 
zt* (E) contains now an infinite product over m 

00 iSp/h—in /2^p 

^T<E)= n na-^^nT^)- (142) 

Tp<T*m=l M^Pl 

As above we shall assume that all periods of primitive periodic orbits are non- 
commensurable. Therefore e'^'^^'^ can be considered as independent random 
variables but one cannot ignore the existence of products over m in Eq. ( |142| ) 

z*T,{E + ei)zT4E + e2) yr . p w ^4-^1 

< z,,{E + e,)z*,4E + e2) < R.i<f>.) >> (143) 



where 



and 



and 



<Rp>= Rp{<Pp)^, (144) 

Jo 271 

RM = n ^'"^ ^, (145) 

iii (1 - A,,^e^<^^W^')(l - A,,^e-<^^-i'') 



Ap,^ = |Ap|-i/2A;-, 0p = Sp(E)/;i - 7r/ip/2, r« = T^e,. (146) 

Though it is possible to calculate the integral < Rp > hj the residues 

as has been done in the previous Section it is more convenient to use the 
so-called g-binomial theorem (p2| p. 7) 

£ ^^z- = (147) 

Here 

(a; g)„ = (1 - a)(l - ag)(l - aq^) ... (1 - ag^'^), (148) 
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and 

Denoting 
one obtains 

Similarly 



a 



;?)oo= n(i-0- (149) 



n=0 



(2) 

^ (1 - A,,^e-^^+'-^ ) ^ (g^^; g)oo 



where z' = |Ap|"^/^e*'^f+*^p . 

Using the g-binomial theorem these expressions can be represented as 
follows 

„ [i Ap.^e _) ^ ^ (153) 

and 

Now < -Rp > is the integral of the product of these two expressions and in 
the diagonal approximation only terms with the same power will contribute 
and finally 

<Rp>=j: V^i/", (155) 

where y = \Ap\ '-e'^^'^p This sum may be expressed through the so- 



called g-hypergeometric function ||22 

< Rp >=2 0l(e-*^''^e-*^''^ A;^ A;\ \Ap\-'e'^-'), (156) 
where we take into account that r^^^ — rj^"^^ = Tp{e). 
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The total mean value equals the product over all periodic orbits till T* . 
The divergent part comes only from n = 1 term 



n (1+ 



Tp<T* 



2)). 



(157) 



But the divergence of this product coincides with the divergent part of the 
expression 



n 



Zp{ie) 



where 



The product 



ZJs) = 1 - 



\Ar 



U <Rp> 

Tp<T* 

converges when T* — >■ oo and 



Zp{ie) 



R2{e) 



47r2 



-|7-'^d(^e)|^$(°^^)(e)+c.c., 



where 



<|.(°^'=)(e) = [] < > 



Zp{ie) 



Zp{0) 



(158) 
(159) 
(160) 

(161) 
(162) 



and Zci{s) is a classical zeta function (see the footnote at page |TU]) defined 
when Res < by the product over all ppo 



Zri(s) 



n(i 

ppo 



Let us fix the maximal period T* from the condition 

1 



d n ^p(o) 



Tp<T* 



27l\ 



(163) 



(164) 



where 7 is the residue of Zci{s) at s = 1 (which existence is a consequence of 
Hannay-Ozorio de Almeida sum rule) 



7 = lim sZci{s). 

s—>0 



(165) 
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As above this renormalisation defines T* of the order of Th and ensures that 
when e — >^ -R2(e) tends to the oscillatory part of the GUE result (|37|) 

R^^^^\e) = -l—{e^^'^' + e""). (166) 

In principle any zeta function whose divergent part at real e coincides with 
Zci{ie) can be used in the above expressions. 



5 Random matrix universality 

This section is devoted to another method of semi-classical calculation of 
off-diagonal part of correlation function based on a different idea. 

In is well known |^ that the standard ensembles of random matrices 
correspond to the following choice of measure in the ensemble 

P(M) = exp(-Tr\/(M)), (167) 

where V{M) is an arbitrary function. The important universality in the 
random matrices ensembles is the fact that the unfolded distribution does not 
depend on the explicit form of this function |^ provided that it corresponds 



to the so-called definite momentum problem 

For unitary ensembles all correlation functions have the form [j^ 

Rn{Xi, . . . ,Xn) = det\KN{Xi,Xj)\ij=i^,„^n, (168) 

where the kernel 

N 

K^ix, y) = e-^W/2-v^fe)/2 ^ p^{x)p^{y), (169) 

n=l 

and Pn{x) are polynomials of degree n — 1 orthogonal with respect to the 
measure exp(— V"(x)) 

Pn{x)Prn{x)e~^^'''^dx = 5nm- (170) 

The 'semi-classical' asymptotic of orthogonal polynomials leads to the fol- 
lowing expression for the kernel K{x,y) = K]sr{x,y) when N ^ oo 

7r(x - y) 
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where N{x) = p{x')dx' is a mean staircase function related to V{x) by the 
Dyson equation 

P^^^ dt = \v\x). (172) 



x-t 2 

The mean density p(x) does depend on the form of V{x) but if 

e = X — y characteristic scale of the potential V{x) 

one can expand the difference N{x) — N{y) into powers of e, N{x) — N{y) = 
de, and all correlation functions will be functions of de which ensures their 
universality after rescaling. 

Now we shall use this fact in addition to the trace formula. 

Let us assume that we know all periodic orbits up to period T* but, 
contrary to what was used in the previous sections, this maximal period will 
be much smaller than the Heisenberg time but still much bigger than pure 
classical time 

U < T* < tH. (173) 

Below we assume that the maximal period is of the order of the Ehrenfest 
time 

T* ^ tE. (174) 
According to the trace formula the density of states has the form 

d{E) =d{E)+r]{E), (175) 

where as above 

d{E) =d+ H ^^A^^^e''^s,in~.^,i2) ^ ^yj^^ 

Tj,<T* n=l ^ 

and ri{E) is (unknown) part of the density constructed form high-period 
orbits. 

Let us now try to construct a random matrix ensemble which has the 
mean density exactly equals d{E). In principle the necessary potential can 
be computed from the Dyson equation ( [L72| ). But the explicit form of this 
potential is irrelevant as all correlation functions depend only on the mean 
staircase function N{x) = d{x')dx' which is supposed to be known. 
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Hence, the 2-point correlation function in such ensemble takes the form 

~ , , ~, , ~, , sin'^ ( 7i( N( E + x) - N(E + y))) , , 
i?2 X, y) = d{E + x)d{E + y) ^ ^ \ , 

d{E) in this formula is the mean (over constructed ensemble) density of 
eigenvalues and it is not a constant but includes periodic orbits with period 
ut to T* and consequently it has oscillations in the scale smaller or equal 
T*. To be consistent with the standard definition of the 2-point correlation 
function (|38D it is necessary to perform an additional smoothing over over an 
energy window, AE, much bigger than h/T* 

—AE > 1, (178) 

Denoting this average by the brackets < ... >, the 2-point correlation func- 
tion 

R2{x,y)=<d{E + x)d{E + y) ^ ^ _ ^ —>■ (179) 

vr [X y ) 

Separating the smooth and oscillatory parts one gets 
Rf"'^ {x, y) = < d{E + x)d{E + y)> 



27r2(x — yy 



where 

1 



A(x, y) = . ^ exp(-47r' < N[E + x)N[E + y) >), (181) 
[x - yy 



and 



47r^(x — y)^ 

As above we assume that all periodic orbit actions in the explicit expression 
for d{E) are non-commensurable. Because we chose T* to be of the order of 
the Ehrenfest time, one can perform the smoothing over energy window in the 
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strict diagonal approximation. But all the expressions which have to smooth 
are exactly the same as were above. The only difference is that the limiting 
period T* is of the order not of the Heisenberg time but of the Ehrenfest time 
and there is a factor (x — Using the same transformations as above one 
concludes that 

A(e) = 4l^d(^e,r*)|2$('^-3)(e), (183) 

and 

iZ,zfeT*) (184) 



47r2e2 



Zci{0,T* 

Here e = x — y and functions <|)('^*'*5)(e) and ^^°^-^\e) are defined exactly as 
Eqs. ( |50D and ( |162| ) except that they are truncated at T*. Because these 
products converge and the value of T* is assumed to obeys inequalities ( |173| ) 
one can go to the limit T* —>■ oo and the resulting expressions will coincide 
with Eqs. (^) and ([1621) . 

The function Zci{s,T*) is the truncated product 

sT 

Zciis,n= n (1-^)"'- (185) 

Tp<T* M^pI 

Because T* is such that this product includes a large number of terms one 
can consider the behavior of the product in the limit T* — > oo. 

Let us perform the following formal transformations (which can be done 
rigorously by adding a small imaginary part to e 

Zci{ie,T*)_l 1 ~ pfa ^7{-\ nQa\ 
11 ; 7^ = -Zci{ie), (186) 



where Zci{ie) is the classical zeta function ( |163| ) defined as the product over 
ppo with arbitrary period and 

7=^ n n -4^- (187) 

Tp<T* |Ap| Tp>T* i - ^ 

But 

1 



n ^ — ^=iim n 



Tp<T* ^ |Ap| Tp<T* 1 - 

isT 

= limZ,,(^.)( n (1-^)- (188) 
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Therefore 



When T* ^ oo 



(1 _ SlllL) 

7 = ehmZ,K^s) n i^l - (189) 



Tp>T* MVI Tp>T* MVI 



Using Hannay-Ozorio de Almeida sum rule one concludes that when 
T* ^ oo this expression is close to 



exp(- /_ —dT). (191) 

Finally for large T* 

7 = elimZ,,(zs)exp(- / -TfrdT + -rdT). (192) 



As the limiting period has been chosen such that eT* ^ 1 (and of course 
sT* <^ 1) one can use the approximation 



roc g^"^-* 

/ —dT ln(T*e) + const, (193) 

JT* T 

Therefore 

-f = \imsZ,i{is), (194) 

and as only |7| is important Eq. ( |184| ) coincides with Eq. (|161|) derived by 
completely different method. 



6 Riemann-Siegel form of density of states 

In this Section we propose one more method of computation spectral statis- 
tics based on ideas different from the ones discussed in the previous Sections. 
For simplicity we consider only the Riemann case. 

It is known that the best method of calculation of the Riemann zeta 
function is the using the Riemann-Siegel representation of it [Q, (called 
in simple cases the approximative functional equation). 
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Heuristically it can be derived as follows. Let us divide the sum over 
integers in the definition of the Riemann zeta function into two contributions 



oo -I TV -I oo -1 

C(l/2-.E) = 5:^e^^^'^" = 5:^e^^"^"+ Y: ^e^^''^^"- (195) 

n=l n=l n=N+l 

In the first sum the summation is done over integers up to a (large) integer 
and the second sum includes all terms with n > N. Using the Poisson 
summation formula one can rewrite the latter in the form 



oo 1 oo „oo 1 

y — -e'^^""= y / — -e-2~+'^''^"rfn, (196) 



and it is naturally to apply the stationary phase method to compute the inte- 
gral. The dominant contribution comes from vicinities of points Ugp obeying 
the saddle-point equation 

2Tim = — , (197) 

or 

= (198) 

but this contribution exists only if this stationary point lies between and 
oo which is equivalent to the following restriction on m 

0<m< N*, (199) 

where 

N* = (200) 

27iN ^ ' 

Computing the second derivative of the phase and performing the integration 
in the Gaussian approximation one gets 

oo -| _ A''* -| 

y J_e^si-- = e2-^(^) V-!-e-^^"^", (201) 

where N{E) is the mean staircase function given by Eq. ([T7|). (On more 
careful analysis see [pH]). 
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The same result can be obtained without calculation (but less rigorously) 
by using the functional equation of the Riemann zeta function. It is known 
that this function obeys the relation ^ 

C(l/2 - lE) = e2"*^(^)C(l/2 + iE), (202) 

which is a necessary condition that the complex valued function ({s) has 
zeros on a line. Substituting the formal expansion 

C(l/2±zE) = ^^e±^^i"" (203) 

n=l ^ 

into this relation and comparing terms in two parts one concludes that 

E J^e-^^'^ = e-^(-) E r^e---", (204) 

lnn>Ti lnn<T2 "' 

provided Ti + T2 = Th where the Heisenberg time Tjj = 2TTd{E). One can 
also get a slightly more general relation 

N2 1 E/{2nNi) . 

V - p27ri7V(£;) ^ (20^) 

2-^ r,l/2-iE ~ ^ r,in+iE- ^'^^^> 

n=Ni n=E/i2TrN2) 

This type of bootstrap where a sum of large integers is proportional to a sum 
of the short ones is the cornerstone of the Riemann- Siegel method and it is 
very useful in calculation of zeta functions. Our purpose will be to use it to 
obtain an information about sums of primes. 

Let us split the density of Riemann zeros (|T^) into two contributions 

S'"''\E) =d,{E) + d2{E). (206) 

The first one includes the sum over all prime up to a certain integer N 

1 °° 1 

diiE) = E E w^e^^""^^ + CO. (207) 

p<Ar „=i P 

and the second term contains all primes bigger than N 

1 °° 1 

ME) = E E w^e^''""^^ + c-c- (208) 

p>N n=l P 
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The main ingredient of following discussion is the sieve representation of 
prime numbers. Namely we shall use the following relation 



fip) = H /(^) - fiPj^) + Y fiPjPku) + ... 

P>N n>N pj,n>N/pj P3 ,Pk\n>N / {pjpk) 

= E E fink)^i{k), (209) 

k n>N/k 

where ii{k) is the Euler function. 

This relation is just the manifestation of the fact that when one substi- 
tute instead of a sum over prime numbers the sum over integers one has 
first to subtract from this sum all numbers which are proportional to short 
primes but all numbers which are proportional to product of two primes are 
subtracted twice, therefore, it is necessary to add integers proportional to 
the product of two primes etc. This type of inclusion-exclusion principle is 
exact provided one is restricted to a finite collection of primes. 

Applying this formula to f{n) = l/n^/'^~^^ one gets 



E ^1/2-iE E E ^1/2-iE l.l/2-iE' 

(210) 

p>N P k=l n>N/k ^ ^ 



Using the approximate functional equation ( p04| ) for the sums over n in this 
formula one gets 



N2 1 00 Ek/{2TTNi) . /,x 

_ 2niN(E) A'-l'^J 

p=Ni P^' k=l n=Ek/{2nN2) 



^ ^1/2-iE ^1/2+iE f^l/2-iE' K^^-^J 



Taking into account that the total 'period' of an individual term in the double 
sum is t = Ink — Inn this expression can be transformed into the following 
form 



^ „2mN{E) I ^ 1 ^l{k) 

P=Ni P ' \k=\n=\ "' ' ' / TH-\nN2<t<TH-\nNi 



^1/2-iE ^ \^ ^ j^l/2+iE J^l/2-iE I ' V^-^^^ 



where the symbol ti < t < t2 means that one has to expand the expression 
in the brackets into the Fourier series on exp{—iEt) and to take into account 
only terms with period ti < t < t2. It is important to note that if InA^i is 
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of the order of Th the summation in the above formula includes only small 
number of terms. But 

Y ^ I (213) 

and these sum can be rewritten as follows 

V ^_ - e^^^^(^) f ^ilZ^±M^ (214) 
2^ r,l/2-iE-^ \c(l/2-iE) ' ^ ^ 

Ni<p<N2^ \SV-^/^ ''^J / TH-hiN2<t<TH~hiNi 

Using the prime number theorem the main part of the sum of large primes 
in the trace formula for the density of the Riemann zero can be rewritten in 
the form 

p=Ni Z^Ni<p<N2 ^ Ni<p<N2 

Let us use the inclusion-exclusion principle for the denominator 

E 1 = (AT, - AT,) (1 - E 1 + + . . .) = (iV2 - A^i) - ^ (216) 

p=Ni P ^ Pi,P2 ^1^2 p P 



Therefore 



0<t<Tsr-lniV 



and the density of Riemann zeros takes the form 
1 °° 1 

diE)- -^EWEw^^'^'^^'^^ (218) 

p<N n=lP 



\P \ ^ '\ f ^/0<t<TH-lnN 



+ C.C. 



This formula can be considered as a Ricmann-Siegel look-like ressumation of 
the trace formula. It expresses the (unknown) sum over large primes through 
small ones. 
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In particular this formula is useful for the computation of statistics of 
Riemann zeros. Direct application of the definition (pSj) gives 



i?2(e) = Rt"'^ (219) 

^2mdt (^-y _ p-l/2+i{E+ei)^^(^-^ _ p-l/2-i{E+e2)-^ 

47r2 ^ n _ p-l/2-i{E+e^)^^l _ p-l/2+i(£;+£2) ) ( 1 _ p-l)2 ^ " 



Here r'^^"'^^ is the same as in Eq. (^) except that the summation includes 
primes up to A^. Because N is assumed to be large and the sum in (^) con- 
verges conditionally at real e one can go to the limit N oo and Eqs. ([M])- 
( |76D will be recovered. The second term, corresponding to R^^'^\ comes from 
the summation over large primes and it has the form as in Eq. (|130|) (but 
with correct pre-factor). Repeating the same calculation as in Section ^TT] 
and taking into account that the final product now converges one will get 
Eq. (|136|) with the correct value ( |139| ) of the pre-factor. 

7 Conclusion 

The above-discussed methods demonstrate how, in principle, the existence of 
the trace formula and certain natural conjectures about the distribution of 
periodic orbits (or primes) combine together to produce universal local statis- 
tics. For systems without the time-reversal invariance the assumption that 
periods for all (or at least for most) periodic orbits are non-commensurable 
leads to the GUE statistics (for 2-point correlation function). 

In all cases the main formulas may be written in the following forms 

i?2(e) = i?f '^3)(e) + R^^"\e), (220) 
where the diagonal part 

Rf^'\e) = -^^ln(|^c/(6)P<f('^^"^ne)), (221) 
and the off-diagonal part 

R2"\^) = -\e2-'^"^|-Z,Ke)P$(°^^He). (222) 
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Here Zci{t) is a classical zeta function and <|)('^*'^5)(e) and ^^"^f^e) are certain 
convergent products over periodic orbits. 
For dynamical systems 

= 11(1 - exp(zerp)/|Ap|)-\ (223) 
p 

For the Riemann case 

Z,l{e) = C{l-^E). (224) 

In all cases Zci{e) ^ 7/e when e ^ and $*^°-^'^)(0) = 1 which ensures the 
GUE value for small e. 

The close relation between diagonal and off-diagonal terms (first observed 
for disordered systems by Andreev and Altshuler |^B[) suggests the existence 
of a certain unified principle. The best candidate for it is the 'unitarity' 
property of the trace formula, namely, that the distribution of periodic orbits 
should be such that the corresponding eigenvalues will be real. In some 
sense certain long-period orbits are connected to the short ones and the 
investigation of this connection may clarify the origin of universal spectral 
statistics. 

The interesting question remains what conjectures about periodic orbits 
are necessary to obtain correlation functions for systems with time-reversal 
invariance where almost all periodic orbits appear in pairs with exactly the 
same action. 
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